Fermiology of Cuprates from First Principles: 
Prom Small Pockets to the Luttinger Fermi surface 

L. Hozoi, M. S. Laad and P. Fulde 

Max- Planck- Institut fiir Physik komplexer Systeme, Nothnitzer Str. 38, 01187 Dresden, Germany 

(Dated: February 2, 2008) 



Fermiology, the shape and size of the Fermi sur- 
face, underpins the low-temperature physical properties 
of a metal. Recent investigations of the Fermi surface 
of high-Tc superconductors, however, show a most un- 
usual behavior: upon addition of carriers, "Fermi" pock- 
ets appear around nodal (hole doping) and antinodal 
(electron doping) regions of the Brillouin zone in the 
"pseudogap" state. With progressive doping, 6, these 
evolve into well-defined Fermi surfaces around optimal 
doping (Sopt), with no pseudogap. Correspondingly, var- 
ious physical responses, including d-wave superconduc- 
tivity, evolve from highly anomalous, up to Sopt, to more 
conventional beyond. Describing this evolution holds the 
key to understanding high-temperature superconductiv- 
ity. Here, we present ab initio quantum chemical results 
for cuprates, providing a quantitative description of the 
evolution of the Fermi surface with S. Our results con- 
stitute an ab initio justification for several, hitherto pro- 
posed semiphenomenological theories, offering an unified 
basis for understanding of various, unusual physical re- 
sponses of doped cuprates. 

Understanding high-temperature superconductivity in 
quasi two-dimensional (2D), doped copper oxides remains 
one of the most challenging problems in condensed mat- 
ter physics. In spite of varying structural and chemical 
details, the phase diagram of the high-Tc superconduc- 
tors (HTS's) is seemingly remarkably universal: the un- 
doped compounds with nominally one hole per Cu site 
are Mott insulators (MPs) due to strong electron-electron 
interactions [ij. Upon addition of charge carriers (dop- 
ing), the cuprates turn into d-wave superconductors [d- 
SC) at low temperatures, T <Tc iJj. 

The "normal" state for T>Tc is actually very abnor- 
mal, and radically undermines the conventional Landau 
theory of Fermi liquids (PL's). In the so-called under- 
doped (UD) regime, (5<Cl, a d-wave pseudogap (c?-PG) 
characterises the normal state (H- Whether this d-PG 
state is the precursor of d-SC at lower T or its competi- 
tor is hotly debated 0, [1]. Around optimal doping, a 
"strange metal" phase, with most unusual singular re- 
sponses Q is clearly revealed: this is the celebrated non- 
FL metallic state that has been investigated for twenty 
years In the overdoped (OD) regime {5>6opt), low-T 
FL behavior seems to be smoothly recovered. 

Very recently, notable improvements in sample quality 
as well as measuring technology have finally allowed ac- 
curate mapping of the actual dispersion of the quasiparti- 



cles (QP's) and the Fermi surface (FS) of HTS's. Specif- 
ically, angle-resolved photoemission (ARPES) JE H, 0] 
and quantum oscillation (SdH) techniques [1, [J reveal 
crucial, hitherto unmapped features of the evolution of 
the (renormalised) QP dispersion as a function of dop- 
ing. Thus, these works open up the possibility, for the 
first time, of unearthing the link between the electronic 
structure and physical responses of cuprates in micro- 
scopic detail as a function of S. Both ARPES and SdH 
measurements reveal a full FS consistent with conven- 
tional band-structure calculations for 6 > Sopt @, fiot - 
However, in the UD regime, the small "Fermi" pockets in- 
ferred by SdH experiments are in deep conflict with Lut- 
tinger's theorem. For hole doped samples with (5 = 0.1, 
the SdH results yield a carrier concentration xsdH = 0.15 
Q. Under the same conditions, the low-T Hall constant 
is electron-\ik.e . How can this come about ? Existing 
theoretical calculations cannot resolve this issue satisfac- 
torily. And yet, this finding points toward a glaring dis- 
crepancy in our understanding of the electronic structure 
of cuprates. In light of these findings, a consistent the- 
oretical scenario aiming to describe the unique physics 
of HTS's must now base itself upon the appropriate, col- 
lective excitations stemming from the observed, detailed 
shape and size of the FS. 

Here, we study the dispersion of the lowest hole 
and electron-addition states, as well as the evolution of 
the renormalised FS with doping. Using an ab initio 
wavefunction-based formalism, we describe these with 
quantitative accuracy vis-a-vis recent ARPES and SdH 
measurements. Implications of our findings for other ex- 
periments, as well as their connection to earlier model- 
based and semiphenomenological theories, are discussed 
in detail. Our findings lend credence to the view 
that the unique properties of cuprates are those of a 2D, 
doped MI. 

THEORETICAL FRAMEWORK 

The correlation-induced renormalisation effects on the 
valence and conduction energy bands are remarkably 
strong in cuprates. Early attempts to describe these ef- 
fects were based on the t-J model [ij and indicated the 
crucial role played by the strong antiferromagnetic (AF) 
couplings in reducing the effective bandwidths. In the 
three-band context, it was suggested that a doped oxy- 
gen hole would induce short-range ferromagnetic (FM) 
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correlations between adjacent Cu sites [l3|. If the extra 
hole delocalizes over all four equivalent ligands of a given 
CUO4 plaquette Q , these FM correlations would involve 
Cu sites on five plaquettes as shown in Fig. 1. Since 
the mobility of this entity is expected to be small, it is 
often referred to as a FM "spin polaron" . 

An accurate investigation of the structure of such com- 
posite objects calls for methods that allow an unbiased 
treatment of the various (competing) interactions in the 
Cu02 plane. A fundamental point underlying the physics 
of cuprates is the interplay between electron localization 
effects as a result of strong repulsive interactions and 
band-like behaviour as a result of translational symme- 
try and inter-site orbital overlap. Standard band theories 
based upon the density-functional model (and the local 
density approximation, LDA) mainly emphasize the lat- 
ter aspect. Though LDA provides rather good results for 
weakly correlated solids, its limited ability to describe 
correlated d (and /) electrons is well-documented. With 
the advent of dynamical mean-field theory (DMFT), this 
basic conflict has been partially resolved JJiJ. In partic- 
ular, much progress in describing the FS's of real materi- 
als, along with their one-particle spectral functions, has 
been possible. However, this is still some distance from 
being a totally ab initio approach, since the actual corre- 
lations are approximated by local (Hubbard) parameters. 
Use of constrained LDA to estimate these parameters 
entails an uncertainty of the order of 20%, while their 
self-consistent estimation within LDA-t-DMFT is fraught 
with insurmountable problems [l^ . 

An alternative approach bases itself on state-of-the-art 
quantum chemical (QC) methods flS\. In molecular sys- 
tems, wavefunction-based quantum chemistry provides a 
rigorous theoretical framework for addressing the elec- 
tron correlation problem fl9'|. A real-space, QC-based 
treatment is then a natural starting point in dealing with 
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FIG. 1: Ferromagnetic correlations (green) among Cu sites 
around a doped oxygen hole (violet) as inferred from CASSCF 
calculations. Farther spin couplings are antiferromagnetic 
(red). The motion of the oxygen hole is coherent when the 
FM spin polarization "cloud" at adjacent Cu sites moves sol- 
idarity with the hole. The O ligands are at the vertices of the 
square plaquettes, the Cu ions are shown as grey dots. 



Mott physics in d-metal solid state compounds. As shown 
below, the k-dependent energy bands can be recovered at 
a later stage after rigorously accounting for the ubiqui- 
tous strong short-range correlation effects. 

The strategy is to use a sufficiently large cluster, C, 
cut out from the infinite solid and properly embedded 
in some effective lattice potential, capable of describing 
these crucial short-range correlations accurately. The 
presence of partially filled d electron shells requires a 
multiconfiguration representation of the many-electron 
wavefunction. The complete-active-space (CAS) self- 
consistent-field (SCF), CASSCF, method ^ provides 
precisely such a framework (see Methods for details). It 
can, for example, describe spin correlation effects such as 
the Anderson superexchange in Mi's 12111 and the double- 
exchange in mixed- valence systems [22]. For undoped 
cuprates, with formally one 3dx2-y2 electron per Cu site, 
the CAS wavefunction is similar to the variational wave- 
function used in numerical studies of the 2D, one-band 
Hubbard model However, all one- and two-particle 
integrals are computed here in a totally ab initio way. 

The dispersion of d-like states on a square lattice is 
given by the following relation: 

e(k) = — 2t (cos k^a + cos kyo) + At' cos k^a cos kyG 
— 2t"(cos 2kxa + cos 2kya) , 

where t, t' , t" are the hopping integrals between nearest- 
neighbor (NN), second-NN and third-NN sites and the 
effective site is one CUO4 plaquette IJ]. LDA calcula- 
tions yield t values of 0.4-0.5 eV and a ratio between the 
NN and second-NN hoppings t'/<«0.15 for La2Cu04 and 
0.33 for TlaBaaCuOe Q. In contrast, the CASSCF cal- 
culations predict a renormalised NN hopping t = 0.135 
eV in La2Cu04 [1^1 ■ Here, we describe how the de- 
tailed QP dispersion can be obtained with quantita- 
tive accuracy, for both hole and electron-addition states. 
The effective hoppings are computed by using the over- 
lap, Sij, and Hamiltonian, iJjj, matrix elements between 
(iV=pl)-particle wavefunctions having the additional par- 
ticle (hole or electron) located on different plaquettes 
(«,j,...) of a given cluster. Each of these (A^^pl) wave- 
functions, \^^'^^), is obtained by separate CASSCF op- 
timizations. A similar scheme was previously applied to 
simpler, noncontroversial systems such as diamoiid, sili- 
con and MgO 25, 2^. It accounts for both charge [25l[26| 
and spin polarization and relaxation effects in the nearby 
surroundings, see Fig. 1. For degenerate (i.e., Hii = Hjj) 
(TVtI) states, t = (e^ - e,)/2 = {H,j - S^jHu)/{l - Sf^), 
where and ej are the eigenvalues of the 2x2 secu- 
lar problem. For non-degenerate {Ha ^ Hjj) states, 
t = l/2[(e, - e,)2 - {H,, - Hun/' 
H, 

method [27|. Correlation effects beyond CASSCF on the 
onsite matrix elements Ha are calculated by multicon- 
figurational second-order perturbation theory, CASPT2 



The Sij and 



tj terms are computed using the State-Interaction (SI) 
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[28|. The short-range magnetic correlations are included 
in our clusters by adding extra CUO4 units around those 
plaquettes directly involved in the hopping process, see 
Fig. 2. 



QUASIPARTICLE BANDS AND ARPES 

We investigated both, the p-type HTS La2Cu04 and 
the n-type HTS SrCu02- Renormalised hopping matrix 
elements (ME's) for Zhang-Rice (ZR) ^ type states in 
La2Cu04 and electron-addition states in SrCu02, involv- 
ing neighbors up to the third order, are listed in Table 
I. For comparison, "unrenormalised" (or bare) hoppings 
were also computed, by imposing a FM arrangement of 
spins at the nearby Cu sites, i. e., a FM lattice. The bare 
hoppings are substantially smaller for the (N+l) Cu d^'^ 
states because the Cu 3d functions are more compact as 
compared to the O 2p orbitals. 

In effective one-band models [H, , the ZR p-d state 
is regarded as a vacant, or unoccupied, d-like site. Con- 
sequently, there is no renormalisation of the second-NN 
and third-NN hoppings t' and t" because these connect 
sites of the same magnetic sublattice. In contrast, the in- 
terplay between short-range FM correlations and longer- 
range AF couplings (see the sketch in Fig. 1) produces 
large renormalisation effects for all hopping ME's in our 
approach. For t' and t", in particular, the hopping of 
the 2p hole implies coupled, Cu and O spin "flips" on 
the two plaquettes directly involved in the hopping pro- 
cess. Nevertheless, spin correlations decay rapidly with 
distance and so the renormalisation effects are less dras- 
tic for the third-NN ME, t" . We thus find i" ~t/2 and 
t' <^t, a rather remarkable result. 

For the description of the electron-addition c?^" states, 
an effective one-band model is seen to be justified. As 




shown in Table I, in this case only t is substantially 
affected by nonlocal spin correlations. Hence, for the 
(iV-l-1) states, the renormalized hoppings satisfy t' ~t, 
an equally remarkable result. The particle-hole asymme- 
try, see Fig. 3, is now readily understood from the very 
different t' /t and t" /t values for the (A^-l) (ZR-like) and 
(N+l) (Cu d^°) bands. 

Doping of the Cu02 planes is achieved by chemical sub- 
stitution in the "reservoir" layers. The dopant carriers 
must quantum mechanically tunnel from the reservoir to 
the planes: this necessarily involves the apical O 2pz, Cu 
3(ij,2 and Cu 4s orbitals, causing additional renormalisa- 
tion of the planar QP's from these apical charge-transfer 
interactions. Hence, we extended our calculations to in- 
clude configurations where an electron is removed/added 
from/to the Cu 3dz2 or Cu 4s orbital. Only the Cu 3^22 
(TV— 1) state gave rise to considerable renormalisation 
of the planar QP dispersion. While the onsite mixing 
between the ZR and [N—l] configurations and NN 
hopping between degenerate dz2 hole states are negligi- 
ble, the inter-site off-diagonal hopping is large. With 
sets of orbitals individually optimized for each (A^ — 1) 
state, this off-diagonal ME is = 0.20 eV, larger than 
the value reported in ref. [sB], where the lowest d^z hole 
state was expressed in terms of orbitals optimized over 
an average of several excited states involving different 
couplings among the nearby Cu spins. In k-space, the 
hybridisation ME between the ZR and d^^ bands reads 
7,„(k) = tm{coa k^a — cos ky a). Further, onsite, the ZR 
and dz2 (N—l) states are separated by an energy Ae. It 
turns out that the correlation-induced corrections to the 
CASSCF energy separation are substantial, changing this 
quantity from 0.60 eV ^ to Ae = 1.70 at the CASPT2 
level. Such corrections are usually small, for both the 
onsite relative energies [m and hoppings [Ulii]. Nev- 
ertheless, corrections as large as 0.9 eV have been found 
before for the relative energy of the state of the d^ 
manifold in NiO j2ij, for example. 

It is now straightforward to diagonalise the k- 
dependent 2x2 matrix. 



ezR(k) 
7m (k) 



7m (k) 
2(k) + Ae 



FIG. 2: Sketch of the finite clusters employed for the cal- 
culation of the effective valence-band and conduction-band 
hoppings. 



to yield the renormalised bands. This constitutes a non- 
trivial extension of the three-band Hubbard model, where 
the additional renormalisation from the apical link is not 
considered. The resulting dispersion of the ZR-like band 
is plotted in Fig. 3 and shows excellent agreement with 
the dispersion of the lowest ARPES band reported for 
La2Cu04 by Ino et al. [3l|. In particular, the flat dis- 
persion around (0,7r), the maximum near (7r/2,7r/2) and 
a renormalized bandwidth of nearly 1 eV are all faith- 
fully reproduced in the theoretical results. We have not 
attempted to describe the "waterfall" -like structures ob- 
served recently in ARPES. First, their interpretation [3^ 
and causal link to c?-SC are controversial. Theoretically, 
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TABLE I: Hopping ME's for ZR-like states in La2Cu04 and 
electron-addition Cu d^" states in SrCu02. The bare hop- 
pings were computed by imposing high-spin couplings among 
the nearby Cu sites. For the d^" states, each of the hop- 
pings changes by less than 15% from La2Cu04 to SrCu02 
(not shown in the table). On the other hand, no ZR-like so- 
lution was obtained for the (A^ — 1) states in SrCu02, which 
qualitatively confirms the experimental findings: for in-plane 
lattice constants a > 3.87 A (a = 3.925 A in SrCu02), the 
Cu02 planes do not readily accept holes in the O bands [29j. 



Hopping ME's 



Bare 



Rcnormalised 



"ZR" 

t 
t' 
t" 



t 
t' 
t" 



state 



0.540 
0.305 
0.115 



0.290 
0.130 
0.045 



0.135 
0.010 
0.075 



0.115 
0.130 
0.015 



the study of such structures requires an anal ysis of the 
incoherent part of the spectral function [H, 111. This 
challenging exercise is beyond the scope of our present 
work. For the Cu d^'^ states, a lack of detailed data for 
ARPES lineshapes in n-type cuprates precludes a direct 
comparison between theory and experiment. 

Knowledge of the QP dispersion enables us to study 
the evolution of the renormalized FS as a function of 
doping. Assuming a rigid band shift with doping, an as- 
sumption supported by independent experiments 34, 35| . 
we plot the evolution of the FS for both hole and electron 
doped cuprates in Fig. 4. We simulate doping effects 
by a progressive downward shift (hole doping) and up- 
ward shift (electron doping) of the Fermi energy, Ep- 
Once again, our results show a remarkable agreement 
with the experiment: small hole pockets centered around 
the nodal (N) region [k„ = (7r/2a, 7r/2a) ] comprise the 
"FS" in the deeply UD regime [1, ij. Additionally, for 
slightly higher (5, smaller, electron-]Ske pockets centred 
around the corners of the Brillouin zone are also clearly 
resolved, see Fig. 4(a). With further doping, these pro- 
gressively evolve into a large hole-like FS [5| , implying a 
FS reconstruction close to 5 opt- At a critical value 5 = 5c, 
the FS changes from hole-like {5 < 6c) to electron-like 
6 > 6c), in complete accord with results from ARPES 
5|[. Excellent agreement of the FS vis-a-vis experiment 
36| is also obtained for the n-type cuprates: small pock- 
ets centered around ka„ evolve into a hole-like FS with 
progressive electron doping. To our knowledge, ours are 
the first ab initio results capturing such effects: hith- 
erto, these have been (partially) described within effec- 
tive, one- 



lil . 33 . ST*! or three-band [H, [s^ models with 



BROADER IMPLICATIONS 

Our results constitute an ab initio derivation of the 
"hot-spot-cold-spot" phenomenology also seen in 

cluster-DMFT work on the 2D Hubbard model [13]. 
There, the QP scattering rate is strongly k-dependent 
as the FS is traversed. This is mainfest in our com- 
puted FS: the pronounced dispersion of the N-QP's im- 
plies weaker QP scattering around k„, in marked con- 
trast to the antinodal (AN) region, where strong band 
flattening is indicative of strong QP scattering and very 
short QP lifetimes. Further, the renormalized t, t' , t" 
imply intrinsically frustrated hopping: interestingly, the 
importance of frustrated kinetic energy to the high- Tc 
problem has been discussed at length in the resonating 
valence bond (RVB) model of Anderson 12:|. We show 
that both these seemingly disparate features arise from 
the same underlying microscopic mechanism: strongly 
anisotropic renormalisation of carrier motion by strong, 
short-range spin correlations. 

Interestingly, in a Hubbard-type model, large t' (or 
t") open the door to additional exotic phases, like d- 
wave nematic [iH, d-density wave [i^l and valence-bond 
fi^l ordered phases. These have been invoked as possi- 
ble competitors of rf-SC in various semiphenomenological 
contexts. Our work establishes the intimate connection 
between these putative instabilities and short-range spin 
correlations characteristic of a (lightly doped) MI. 

Our findings have remarkable implications for the in- 
terpretation of a host of experiments probing the un- 
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FIG. 3: (a): The ZR-like electron-removal band for La2Cu04 
in the 2D Brillouin zone, without including the interaction 
with the d^^2_^2 hole state (dashed line) and after including 
this interaction (thick red line). For clarity, the do,^2_^2 band 
is not shown in the figure. The zero of energy is the value of 
the onsite Hamiltonian ME of the ZR state, Hf^. (b): QP 
dispersion for the electron-addition Cu d^" state in SrCu02. 
The reference energy is the value of the onsite Hamiltonian 
ME, Hfi . Units of eV are used in both panels. 
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FIG. 4: (a): 2D colour map for the ZR-like QP dispersion in La2Cu04, including the efltect of the ZR-d^2 hole interaction. 
In a rigid-band picture, the constant-energy (CE) contours (black curves) illustrate the evolution of the FS with hole doping, 
from small (hole) pockets in the nodal region for deeply UD samples to a large, hole-like FS at intermediate dopings and an 
electron-like shape in the OD regime [Bl . For a narrow doping interval in the UD regime, the "FS" is defined by eight contours 
(see the dashed curves) , with both hole-like sheets and electron pockets centred at the corners of the Brillouin zone. The latter 
are related to the small dip in the dispersion at the (tTjTt) point, see Fig. 3(a), and have been also inferred from Hall-effect 
measurements [ij]. (b): Colour map for the electron-addition Cu states in SrCu02. The CE contours reproduce the 
experimentally observed evolution of the FS with electron doping, from small electron pockets in the antinodal region for the 
lightly doped regime to a large, hole- like FS at high doping [3, |3^. The same energy scales as in Fig. 3 are used. 



usual "normal" -state physical responses of the cuprates 
as a function of doping. First, the UD p-type cuprates 
are revealed to be N-metals for small 5, as inferred by 
ARPES [i,|3]- All the "nodal metal" phenomenology % 
can now be provided an ah initio justification in light 
of our results. Coherent nodal QP's also naturally ex- 
plain the "good metal" thermal conductivity in the UD 
region 4j]. Additionally, the SdH quantum oscillation 
frequencies f2 [1, Q now correspond to coherent carrier 
orbits in the "pockets". The Onsager-Lifshitz formula, 
ri = <i>o^/(27r)^, with $o = hc/e the flux quantum and 
A = 27r^a;/4a^ the area of the pocket for x carriers in 
k-space, tells us that ^{x)^x. Our finding of additional 
electron-\\kB pockets around (tt, tt) in the hole doped case 
offers a route toward a resolution of one of the central 
controversies surrounding recent SdH experiments, where 
a^SdH = 0.15 for the 5 = 0.10 cuprates [8|. While this is 
irreconcilable with theories having only hole pockets, our 
finding of additional electron-like sheets can reconcile the 
SdH results with the Luttinger sum rule for ^ss 0.1 
Moreover, the Hall constant Rh{x) is now expected to 
track the evolution of the renormalised FS with doping. 
Depending upon the concentrations of hole and electron 
carriers nh and Ue, with Uh > rie, and the mobilities 
and [lili uiay change sign from hole-like to 



electron-like with T, reconciling the SdH 0, j9|] and Hall- 
effect [llj data. 

Finally, what about superconductivity? Given the 
small number, S, of nodal QP's in the UD regime, d-SC 
will result from pairing of these quasicoherent fermionic 
entities. Without going into the nature of the pairing 



mechanism, the above implies that Tc{S)^5 for 5^ 1. So 
the superfluid density at T = will decrease linearly as d 
is reduced, implying that upon underdoping d-SC will be 
progressively destroyed by order parameter phase fluctu- 
ations, which grow as the MI is approached. This implies 
non-BCS (large) values of 2A/kBTc and strong vortex- 
liquid like effects above for UD cuprates. These have 
indeed been invoked in connection with the anomalous, 
"giant" Nernst effect in UD superconductors How- 
ever, once a full FS develops, d-SC would be expected to 
revert back to a more conventional BCS-like variety [ist . 
with identical scaling for the N and AN gaps. 



CONCLUSIONS 

To summarize, we have implemented a first-principles, 
wavefunction-based calculation of correlated hole and 
electron-addition quasiparticle states in layered cuprates. 
In addition to quantitatively describing the dispersion 
of the ZR-like band, our work reproduces the FS evo- 
lution as a function of doping, in remarkable agreement 
with a host of recent ARPES and quantum oscillation 
experiments. Our finding of large longer-range effective 
hoppings implies intrinsically frustrated carrier kinetic 
energy, in agreement with Anderson's RVB ideas [l^ . 
The very different behavior of hole and electron doped 
cuprates is clearly manifested as originating from very 
different quantum chemical and spin correlation "back- 
grounds". Seen from this perspective, the FS "recon- 
struction" with doping ll|, as well as the famed nodal- 
antinodal dichotomy in the UD systems, are both un- 
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derstood in terms of k-space differentiation of QP states 
in the 2D, doped MI. Phenomenologically, the computed 
evolution of the FS with 6 goes hand-in- hand with the ob- 
served evolution of d-SC from a strongly non-BCS, phase 
fluctuation dominated type, to a more conventional BCS 
type with progressive doping, benchmarking the crucial 
relevance of fcrmiology in cuprates. 

— Methods — 

In the CASSCF approach [l^,!!^], the wavefunction is writ- 
ten as a linear combination of configuration state functions 
(CSF's) |m), |*> = EmC-m \m). These CSF's are spin- 
(and symmetry-) adapted combinations of Slater determi- 
nants (SD's), i. e., eigenfunctions of the operators for the pro- 
jected and total spins. In turn, the SD's are constructed from 
a set of real and orthonormal spin orbitals {(j>p{r,a)}, where 
r and a axe the spatial and spin coordinates, respectively. In 
this work, an initial guess for these orbitals is obtained from a 
Hartree-Fock calculation for an hypothetical Cu 3d^'^, O 2p^ 
closed-shell configuration of the Cu and O species. 

In determining the CASSCF wavefunction, the orbitals are 
variationally optimized simultaneously with the coefficients 
of the CSF's. The orbitals employed for expressing the wave- 
function are thus the optimal orbitals for the state at hand 
and do not introduce a bias toward a particular configura- 
tion. Three different sets of orbitals are used in CASSCF: 

(i) the inactive levels, doubly occupied in all configurations, 

(ii) the virtual orbitals, unoccupied in all configurations, and 

(iii) the active orbital set, where no occupancy restrictions 
are imposed. For undoped cuprates, with formal Cu 3dP and 
O 2p® valence states, the active space would include the par- 
tially occupied in-plane Cu 3d^2_y2 orbitals. Such a CAS 
wavefunction is similar to the variational wavefunction used 
in numerical studies of the 2D, one-band Hubbard model [2^ . 
The main difference is that all integrals, including inter-site 
Coulomb and exchange terms, are computed here in a totally 
ab initio way. In particular, the lower, completely filled lev- 
els, e.g., the O 2s and 2p orbitals, do affect (i.e., screen) the 
actual interactions among the active electrons by readjusting 
themselves to fluctuations within the active orbital space. 

If extra holes are created, the active space must be enlarged 
with orbitals from the inactive group. Each additional doped 
hole requires one orbital to be transferred from the inactive 
to the active space. For the lowest electron-removal state, 
for example, the orbital added to the active space turns into 
a ZR-type p-d composite |l4| in the variational calculation, 
localized on a given Cu04 plaquette. With regard to the 
electron-addition conduction-band states, these turn out to 
have Cu 3d^° character and, in a first approximation, an active 
orbital space including only the 3d^2_y2 levels would suffice. 

So-called dynamic correlation effects [l^ for the onsite ma- 
trix elements Ha were computed using second-order pertur- 
bation theory (the CASPT2 method ^). The Cu 3d and 
O 2s, 2p electrons on five plaquettes (i.e., the ZR plaquette 
and the two apical ligands for that plaquette plus the four 
NN plaquettes) were correlated. All calculations were per- 
formed with the QC software molcas .4(^. For the ions of 
the plaquettes directly involved in the hopping process, all- 
electron basis sets (BS's) of triple-zeta quality were applied. 



These were Gaussian-type atomic-natural-orbital BS's from 
the MOLCAS library [4^, with the following contractions [i9l |: 
Cu (2Isl5pl0d)/[5s4p3d] and O (I4s9p)/[4s3p]. The core 
electrons of the remaining ions of each cluster, see Fig. 2, 
were represented by effective core potentials (ECP's), i.e., Cu 
ECP's plus valence double-zeta BS's [13] and O ECP's with 
triple-zeta BS's [i^. To describe the finite charge distribu- 
tion at the sites in the immediate neighborhood of the cluster, 
we model those ions by effective ion potentials [4^. Beyond 
these neighbors, we use large arrays of point charges (PC's) 
that reproduce the Madelung field within the cluster region. 
Apical ligands are explicitly included in our calculations only 
for the "active" plaquettes. Other apex oxygens are repre- 
sented by formal PC's. That the charge populations of the 
Cu 3d and active O 2p orbitals are not sensitive to the size 
and shape of the clusters we use was shown in ref. [Tsj . We 
employed the crystal structure measured by Cava et al. [H^ 
for La2Cu04 and by Smith et al. [29| ] for SrCu02. 
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